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The guiding-center kinetic neoclassical transport code, XGCO, [C.S. Chang et. al, Phys. Plasmas 11 , 2649 
(2004)] is used to compute the heat fluxes and the heat-load width in the outer divertor plates of Alcator 
C-Mod and DIII-D tokamaks. The dependence of the width of heat-load fluxes on neoclassical effects, neutral 
collisions and anomalous transport is investigated using the XGCO code. The XGCO code includes realistic X- 
point geometry, a neutral source model, the effects of collisions, and a diffusion model for anomalous transport. 
It is observed that width of the XGCO neoclassical heat-load is approximately inversely proportional to the 
total plasma current Ip. The scaling of the width of the divertor heat-load with plasma current is examined 
for an Alcator C-Mod discharge and four DIII-D discharges. The scaling of the divertor heat-load width with 
plasma current is found to be weaker in the Alcator C-Mod discharge compared to scaling found in the DIII-D 
discharges. The effect of neutral collisions on the l//p scaling of heat-load width is shown not to be significant. 
Although inclusion of poloidally uniform anomalous transport results in a deviation from the l//p scaling, the 
inclusion of the anomalous transport that is driven by ballooning-type instabilities results in recovering the 
neoclassical l//p scaling. The Bohm or Gyro-Bohm scalings of anomalous transport does not strongly affect 
the dependence of the heat-load width on plasma current. The inclusion of anomalous transport, in general, 
results in widening the width of neoclassical divertor heat-load and enhances the neoclassical heat-load fluxes 
on the divertor plates. Understanding heat transport in the tokamak scrape-off layer plasmas is important 
for strengthening the basis for predicting divertor conditions in ITERj^ 


I. INTRODUCTION 

Insight regarding heat transport in tokamak scrape-off 
layer plasmas is important for strengthening the basis 
for envisaging divertor conditions in existing and future 
tokamak device^^lf^. The XGCO particle in cell 4D gyro- 
kinetic torodially symmetric code with specified back¬ 
ground diffusion coefficients can be used to study orbit 
width effects on the scrape-off layeJ^. In the research 
presented in this paper, the XGCO code is used to inves¬ 
tigate kinetic neoclassical behavior of the heat fluxes on 
the divertor plates in Alcator C-Mod and DIII-D toka¬ 
maks employing realistic divertor geometry. It is impor¬ 
tant for understanding the dependencies that affect the 
divertor heat-load fluxes in planning experiments and in 
developing new models to describe the physics associated 
with the scrape-off layer region. 

The neoclassical divertor heat-load fluxes are com¬ 
puted and changes in the fluxes associated with the intro¬ 
duction of the effects of neutral collisions and anomalous 
transport are examined. In particular, the dependence 
of the divertor heat-load width on the plasma current is 
studied. The anomalous transport in the XGCO neoclas¬ 
sical code is applied in three different ways: i) Within 45 
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degrees from the midplane, ii) uniformly for all poloidal 
angles and hi) using the Bohm and gyro-Bohm type scal¬ 
ings. 

The discharges examined include an Alcator C-Mod 
discharge and a series of four DIII-D discharges that 
represent a plasma current scaiP. The Alcator C-Mod 
discharge 1100212024 was a discharge in an Alcator 
C-Mod/DIII-D similarity campaign. The total plasma 
current in the Alcator C-Mod discharge is 0.9 MA and 
the toroidal magnetic field is 5.4 T. The discharge is 
an enhanced Dc^ (EDA) high-confinement discharge with 
auxiliary heating of approximately 3.9 MW provided by 
ICRH. In the DIII-D discharges, the total plasma current 
is varied from 0.51 MA to 1.50 MA with a fixed toroidal 
magnetic field {Bt ~ 2.IT), plasma triangularity ((5 
0.55), and normalized toroidal beta {/3^ ^ 2.1 — 2.4). 
The total plasma current and auxiliary heating power 
are given in the Table [I| The plasma density at the top 
of the pedestal varies from approximately 4.5 x 10^^m“^ 
in the high plasma density DIII-D discharge, 132016, to 
2.5x 10^^m“^ in the low plasma density DIII-D discharge, 
132018. 

The divertor heat-load width is found to be broader in 
the simulations of discharges with lower plasma current. 
The neoclassical effects are essential to explain the l//p 
scaling that has been observed in tokamak discharges. 
However, the l//p scaling in Alcator C-Mod discharge is 
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DIII-D Discharges 

EFIT time 

msec 

Plasma current 

MA 

Auxiliary heating power 

MW 

132016 

3023 

1.50 

8.12 

132014 

3023 

1.17 

7.36 

132017 

2998 

0.85 

8.50 

132018 

1948 

0.51 

7.10 


TABLE L The DIII-D discharges analyzed in this study and their parameters. 


found to be weaker as compared to the scaling observed 
in the DIII-D discharges. Thi s result is consistent with 
experimental obser vat ion j ^ ^ I ^ ^ i The effects of neutral col¬ 
lisions on the neoclassical heat-load width are compared 
to the neoclassical heat-load width computed without the 
inclusion of neutral collisions. Neutral collisions at the 
edge of the plasma are shown to have a limited effect 
on the divertor heat-load width. It is found that the 
poloidally uniform distribution of anomalous flux can sig¬ 
nificantly alter the 1/Ip scaling. The gyro-Bohm scaling 
of anomalous transport has somewhat weaker effect on 
the divertor heat-load width 1/Ip scaling than does the 
Bohm scaling of anomalous transport. 

The manuscript is organized as follows: In Sec. [TTj fol¬ 
lowing a brief description of XGCO neoclassical code, the 
computation of XGCO heat fluxes on the divertor plates 
is presented. Sec. |II A| contains a description of the neo¬ 
classical prediction for the Alcator C-Mod and DIII-D 
divertor heat fluxes. The effects of neutral collisions on 
the divertor heat-load width prediction are presented in 
Sec. IIB The anomalous transport model in the XGCO 
code is discussed in Sec. nni The effects of anomalous 
transport on the width of heat flux on the divertor are 
presented in Sec. [TV] and the results are summarized in 
Sec. El 


II. COMPUTATION OF XGCO NEOCLASSICAL HEAT 
FLUXES 

In the XGCO neoclassical cod^^^, the gyrokinetic 
plasma ion and electron guiding centers evolve in time 
within a five-dimensional phase space using a realistic 
magnetic equilibrium and limiter geometry. The code 
uses cylindrical coordinates so that the separatrix and 
X-point region can be easily included in the simulation 
domain. In the XGCO code, the motion of charged and 
neutral particle are followed, and the effect of collisions 
are taken into account using a Monte Carlo approach. 
The XGCO code includes models for the source of neu¬ 
trals at the wall and solves the Hamiltonian guiding cen¬ 
ter equations of motion described in Ref. [TsUISl 

A combination of cylindrical and field-line-following 
magnetic coordinates are used in the XGCO code. Mag¬ 
netic coordinates have the advantage of accurately treat¬ 
ing physical phenomena that are aligned with the mag¬ 
netic field lines. Cylindrical coordinates have the advan¬ 
tage of being able to treat the open magnetic surfaces 



FIG. 1. The XGCO neoclassical heat flux to the outer 
divertor plate (Q?ot) versus distance mapped to mid¬ 
plane (rmid) for the Alcator C-Mod discharge 1100212024 
computed using plasma current scaling factors Csp = 
0.33, 1.00, 1.50, and 2.00. 


of the scrape-off-layer and the closed magnetic surfaces 
of the plasma core with equal accuracy and without the 
divergence associated with magnetic coordinates at the 
separatrix. 

The XGCO neoclassical heat fluxes to the divertor 
plates (Qtot) are computed for the Alcator C-Mod dis¬ 
charge and for the four DIII-D discharges. In order to 
test the sensitivity of the heat flux with regard to plasma 
current, a multiplicative scaling factor Cbp is introduced. 
The scaling factor Cbp is an internal numerical multiplier 
introduced in the XGCO code in order to scale the equi¬ 
librium poloidal flux. The toroidal flux is not modified. 
Thus, the amplification factor Cbp can be also considered 
as a scaling factor for the total plasma current. 

In Fig the computed neoclassical heat flux to the 
outer divertor plate Qtot shown for the Alcator C- 
Mod discharge 1100212024, a discharge that was a part of 
Alcator C-Mod/DIII-D similarity campaign. The plasma 
current scaling factors. Cep, ranged from 0.33 to 2.0. The 
value of Qtot found to be large in the region close to 
the separatrix. A lower peak value of Qtot ^ broader 
width is obtained when Cbp = 0.33 than obtained with 
higher values of Cep. It is seen that as the Cbp scaling 
factor is increased from 1.0 to 2.0, the width of the heat 
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FIG. 2. The XGCO neoclassical heat fluxes to the outer divertor plate versus distance mapped to midplane (rmid) for four 
DIII-D discharges. The plasma current scaling factors used are Cup = 0.33, Cup = 1.0, Cup = 1.50 and Csp = 2.0. The results 
shown in panel (a) are for discharge 132016 with Ip = 1.50 MA; in panel (b) for discharge 132014 with Ip = 1.17MA; in panel 
(c) for discharge 132017 with Ip = 0.85 MA; and in panel (d) for discharge 132018 with Ip = 0.51 MA. 


flux to the divertor plate continues to decrease although 
the change is small relative to the change observed in the 
results obtained for the D3D discharges indicating that 
the Alcator C-Mod discharge has a weaker scaling of the 
divertor heat-load width with plasma current. 

In Fig. results for the heat flux, Qtot, versus rmid 
are presented for the four DIII-D discharges considered. 
The four discharges differ primarily in terms of the total 
plasma current, as shown in Table [ij It is observed that 
an increase in the Cpp scaling factor results in a decrease 
in the width of the heat-load on the divertor plate. For 
values of Cpp >1.0 the peak value neoclassical heat flux 
to the divertor decreases as Cpp increases. This is more 
pronounced in the simulation results for the higher cur¬ 
rent DIII-D discharge, 132017, shown in Fig. 2(a). The 
magnitude of the neoclassical heat flux is found to be 
generally smaller for the low plasma current discharges. 
Similar to the result for the Alcator C-Mod discharge, 
the heat flux in the DIII-D discharges is also found to 


be large in the region close to the separatrix. However, 
the magnitude of the neoclassical heat flux close to the 
separatrix in the DIII-D discharges is smaller than the 
corresponding magnitude obtained in the simulation Al¬ 
cator C-Mod discharge (shown in Fig[^. 

A. Sensitivity of neoclassical prediction of divertor 
heat-load width on plasma current 

In Fig. 1^ the computed XGCO neoclassical heat-load 
width (Adiv) is plotted versus the plasma current scaling 
factor (Cbp) for the Alcator C-Mod discharge and for 
the four DIII-D discharges considered. It is found that 
the heat-load width decreases with increasing plasma cur¬ 
rent, that is Adiv oc I~^. In the Alcator C-Mod discharge, 
the dependence of Adiv on Ip is weaker as compared to the 
DIII-D discharges, shown in Fig. IT- For the base case 
(C'bp = 1-0), the neoclassical divertor heat-load width 
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FIG. 3. The XGCO neoclassical heat-load width (Adiv) versus 
discharge and b) for the four DIII-D discharges. 

in the Alcator C-Mod discharge is found to be approxi¬ 
mately 2.3 mm, which is about 20% below the experimen¬ 
tally observed valu^ of approximately 3 mm. As shown 
below in Sec. |IV[ the anomalous effects typically increase 
the divertor heat-load width bringing the simulation and 
experimental results closer to one-another. 

In Fig.[^ it is seen that the change in heat-load width 
with increasing scaling factor, Cep, is not the same for 
the four DIII-D discharges. For Cep > 0.75, the de¬ 
crease in diverter heat-load width with increasing Cep 
for the two DIII-D discharges with higher plasma den¬ 
sity, 132014 and 132016, is less than the decrease for the 
two DIII-D discharges with lower plasma density, 132017 
and 132018. As indicated below, this difference may, in 
part, be due to the omission of the effects of collisions. 


B. Effect of neutral collisions on divertor heat-load width 

In Fig.[^ the XGCO neoclassical heat-load width in the 
presence of neutrals is shown as a function of plasma cur¬ 
rent scaling factor Cep for DIII-D discharge 132016. The 
dotted curve shows the dependence of the width of neo¬ 
classical heat-load when the effects of neutral collisions 
are not include, the solid curve shows the heat load width, 
when the neutral collisionsare included. The dependence 
of the divertor heat-load width is weakly affected by neu¬ 
tral collisions especially at lower plasma currents (lower 
values of Cbp) where neoclassical effects are found to be 
strong. 


III. FINDING ANOMALOUS TRANSPORT PROFILES 
WHICH YIELD XGCO PREDICTED PLASMA PROFILES 
CONSISTENT WITH EXPERIMENTAL PROFILES 

Anomalous transport in the plasma edge region has 
been analyzed previously, and it has been found that par- 


a plasma current scaling factor Cbp for a) the Alcator G-Mod 



FIG. 4. The XGGO neoclassical heat-load width for with and 
without the neutral collisions as functions of plasma current 
scaling factor Csp for DIII-D discharge 132016. 


tide pinche s may have an important role in the pedestal 
regioifi^lTD. Experimental observations and some analyt¬ 
ical models suggest that anomalous transport in scrape- 
off layer may be intermittent and significantly larger than 
in the pedestal region. Studies of the anomalous trans¬ 
port in the plasma edge region are typically based on 
the analysis of experimental data and include a robust 
model for neoclassical transport. Since particle and ther¬ 
mal fluxes obtained from analysis of experimental data 
include contributions from both the anomalous and neo¬ 
classical transport, it is necessary to deduct the neoclas¬ 
sical transport in order to obtain the correct anomalous 
fluxes and effective anomalous diffusivities. 

Although theory-based models for anomalous trans¬ 
port are available in the XGCO code, a goal of the study 
described in this section is to derive the anomalous effec¬ 
tive diffusivities that can reproduce the experimental pro¬ 
files when these diffusivities are utilized in the neoclassi- 



















5 



¥ 


cal kinetic XGCO code. It is found that strong pinches in 
all channels of anomalous transport are necessary to re¬ 
produce the experimental profiles. The diffusivities that 
reproduce the experimental profiles are then used in the 
computation of the divertor heat-load fluxes. 

The XGCO code utilizes an enhanced random walk 
model for the anomalous thermal and particle transport. 
Unlike the standard random walk model that uses the 
same diffusion coefficients for all transport channels, the 
XGCO code distinguishes between the particle and ther¬ 
mal transport. The model is based on the Monte-Carlo 
method for the advection-diffusion equation in the phase 
space for the anomalous transport: 


FIG. 5. A typical effective particle diffusivity versus nor¬ 
malized poloidal flux (?/;) used in the simulations of Alcator 
C-Mod and DIII-D discharges. The particle diffusivity pro¬ 
file is selected to reproduce the experimental particle density 
profile, the solid line curve in Fig. 
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FIG. 6. The XGCO simulation results of the plasma density 
as a function of normalized poloidal flux ('0) in the DIII-D 
discharge 132016 for the first eight ion transit periods. The 
solid line indicates the experimental plasma density profile. 
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where the anomalous flux 
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The coefficients Da and xa are the anomalous particle 
and thermal transport coefficients respectively. These 
coefficients can be obtained from theory-based models or 
prescribed as an input profiles in the XGCO code. In this 
model, the Monte-Carlo particles advance in RZ space 
depends on both anomalous particle and thermal trans¬ 
port coefficients: 
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In XGCO code, three regions of constant diffusivites 
are implemented. The regions of constant diffusivities 
are separated by two narrow transitional regions that 
use tanh-fit. The levels of anomalous transport in the 
three regions, the locations of transitional regions and 
their widths are adjustable parameters. The particle dif¬ 
fusivity profiles as well as the electron and ion thermal 


diffusivity profiles are independent of one another. These 
profiles are adjusted to find a steady state solution that 
reproduces the experimental profiles. The anomalous dif¬ 
fusivity profiles are selected so that the resulting density 
and temperature profiles remains close to the experimen¬ 
tal profiles for at least eight ion transit periods. For ex¬ 
ample, the effective diffusivity shown in Fig. is such 
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FIG. 7. DIII-D heat fluxes profiles, with the anomalous transport localized near the midplane, versus distance mapped to 
midplane (rmid) for four values of plasma current scaling factor Csp. Panel (a) contains the results for the high plasma current 
{Ip = 1.50MA) DIII-D discharges 132016. Panel (b), the low plasma current {Ip = 0.51 MA) DIII-D discharge 132018. 


that when it is combined with the neoclassical transport 
computed with XGCO, the experimental profiles, shown 
in Fig. are reproduced. It has been found that par¬ 
ticle and thermal pinches play an important part in the 
pedestal regions of all discharges analyzed in this study. 


IV. EFFECT OF ANOMALOUS TRANSPORT ON THE 
WIDTH OF DIVERTOR HEAT FLUXES 

The anomalous transport in XGCO neoclassical code 
applied in three different ways: i) Within 45 degrees 
from the midplane, ii) uniformly for all poloidal angles 
iii) Bohm and gyro-Bohm type scaled transport. 

The heat flux profiles, with anomalous transport 
within 45 degrees from the midplane, are shown in Fig.[^ 
for DIII-D discharges 132016 and 132018 for four values 
of the plasma current scaling factor. In these simula¬ 
tions, the implementation of the anomalous transport in 
the region within 45° from the midplane implies that the 
modes which contribute to the anomalous transport at 
the plasma edge are ballooning in nature. The heat fluxes 
profiles shown in Fig.j^are found to be somewhat broader 
in comparison with the heat fluxes obtained with only 
neoclassical transport (see Panel a and Panel d of Figj^. 
Thus, the introduction of anomalous transport results in 
widening the divertor heat-load width. The ballooned 
anomalous loss (outboard midplane), even without an Ip 
dependence, recovers neoclassical plasma current scaling 
proportional to This scaling in the presence of 

anomalous transport implies that the effects associated 
with anomalous transport that is localized around the 
outboard midplane is not significant. 

In Fig. the divertor heat-load width is plotted as a 
function of plasma current scaling factor Cpp for DIII- 
D discharges 132016 and 132018. In carrying out these 
simulations the assumption is made that the anomalous 
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FIG. 8. The divertor heat-load width as a function of 
the plasma current scaling factor Csp for DIII-D discharges 
132016 {Ip = 1.50MA) and 132018 {Ip = 0.51 MA). The sim¬ 
ulations are carried out using the assumption that the anoma¬ 
lous transport is localized near the midplane. Poloidally uni¬ 
form and /p—independent anomalous transport can mask the 
neoclassical scaling. This can be seen by comparing these 
results with the neoclassical divertor heat-load width for the 
discharges 132016 and 132018 plotted in Fig[^ 

transport is localized near the midplane. The divertor 
heat-load width decreases with increasing Cpp even in the 
presence of ballooned anomalous transport. This implies 
that the neoclassical dependence of the divertor heat-load 
width on the total plasma current is preserved. 

The divertor heat-load width is plotted in Fig. [^as a 
function of the current scaling factor for the four DIII-D 
discharges when the simulations are carried out with the 
anomalous transport applied uniformly for all poloidal 
angles. It is found that for all four DIII-D discharges 
the simulated heat-load width no longer decreases with 
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FIG. 9. The heat-load width (Adiv) versus a plasma cur¬ 
rent scaling factor Csp for four DIII-D discharges. The 
anomalous transport applied uniformly for all poloidal angles. 
A poloidally uniform, Ip independent, anomalous transport 
eliminates the neoclassical Ip scaling behavior. 



FIG. 10. The effect of Bohm and gyro-Bohm scalings of 
anomalous transport on the predictions of the heat-load width 
for the DIII-D discharge 132016. The gyro-Bohm and Bohm 
scalings for anomalous transport have been implemented in 
XGGO neoclassical code. The magnetic field dependent Bohm 
and gyro-Bohm scalings of anomalous transport also recovers 
the neoclassical type scaling. 


increasing plasma current. A poloidally uniform, Ip inde¬ 
pendent, anomalous transport eliminates the neoclassical 
Ip scaling behavior. 

The effects of Bohm and gyro-Bohm scalings of anoma¬ 
lous transport on the predictions of the divertor heat¬ 
load width for the high current and high density DIII-D 
discharge 132016 are shown in Fig. The heat-load 
width computed using Bohm and gyro-Bohm scalings 
of anomalous transport is compared with the heat-load 
width when these scalings are not used. It is found that 


the magnetic field dependent Bohm and gyro-Bohm type 
scalings of anomalous transport recover the neoclassical 
type /p“^ scaling. However, both Bohm and gyro-Bohm 
scalings result in narrower heat-load widths at large and 
small values of plasma current as compared to the heat¬ 
load width when these scalings are not applied in the 
simulations. 


V. SUMMARY 

A goal of this study is to improve the understand¬ 
ing of scrape-off layer thermal transport. Understand¬ 
ing the physical effects that contribute to divertor heat¬ 
load fluxes is important in planning experiments, design¬ 
ing future tokamaks, and developing new models for the 
scrape-off layer region. The heat fluxes and heat-load 
width in scrape-off layer are computed using the kinetic 
neoclassical XGCO code. The guiding-center neoclassi¬ 
cal transport code XGCO includes realistic X-point ge¬ 
ometry, a neutral source model, collisions, and diffusion 
models for anomalous transport. 

An Alcator C-Mod discharge, which was a part of Al- 
cator C-Mod/DIII-D similarity campaign, and a series 
of four DIII-D discharges that represent plasma current 
scan are considered. The neoclassical divertor heat-load 
fluxes and a relation for the dependence of the divertor 
heat-load width on the plasma current are obtained in 
the simulations carried out. In examining the depen¬ 
dence of divertor heat-load width on plasma current, ef¬ 
fects of neutral collisions and anomalous transport are 
taken into account. Changes in the neoclassical divertor 
heat-load fluxes associated with the introduction of the 
neutral collisions and anomalous transport are described. 
Anomalous transport in XGCO neoclassical code is ap¬ 
plied in three different ways: i) Within 45 degrees from 
the midplane, ii) uniformly over all poloidal angles and 
iii) with Bohm and gyro-Bohm type scalings. The effec¬ 
tive particle and thermal diffusivities that are used in the 
simulations are such that when they are combined with 
the neoclassical transport computed with XGCO code, 
the experimental profiles are reproduced. The profiles 
for the particle diffusivity as well as for the electron and 
ion thermal diffusivities are independent of one another 
and are adjusted to find a steady state solution so that 
the resulting profiles remains close to the experimental 
profiles for at least eight ion transit periods. 

The simulations of Alcator C-Mod and DIII-D dis¬ 
charges indicate that the variation of neoclassical heat¬ 
load width on the outer divertor plates with plasma cur¬ 
rent is given by I~^ where a is approximately 0.3 for the 
Alcator C-Mod discharge and 0.8 for DIII-D discharges. 
The computed heat-load width in the simulation of Al¬ 
cator C-Mod discharge is roughly consistent with exper¬ 
imental measurement of approximaltely 3 mm for this 
discharg^i^. The computed value of 2.3 mm is obtained 
in a neoclassical limit when the anomalous effects are ex¬ 
cluded. Inclusion the amomalous effects results in a wider 
















width and brings the head-load width prediction closer 
to the experimental value. The weaker dependence on 
/p found for Alcator C-Mod can be attributed to higher 
collisionality in Alcator C-Mod relative to DIII-D. In gen¬ 
eral, the heat flux is found to be large in the area close to 
the separatrix and the heat flux profiles become broader 
for lower values of plasma current. The neutral colli¬ 
sions do not significantly modify the dependence of the 
neoclassical heat-load width on plasma current. The de¬ 
pendence of the neoclassical heat-load width on plasma 
current is modified when the anomalous transport is ap¬ 
plied uniformly for all poloidal angles. However, The 
neoclassical dependence of the divertor heat-load width 
on the plasma current is preserved when the anomalous 
transport is applied in the region within 45° from the 
midplane. The Bohm and gyro-Bohm scalings of anoma¬ 
lous transport in XGCO simulations results in narrower 
divertor heat-load widths that exhibit the depen¬ 
dence. 
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